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There are five sections in this chapter. They are 
intended to be studied consecutively in four study 
sessions. Section 4 requires the use of the 
computer and Computer Book C. 

The pattern of study for each session might be as 
follows. 

Study session 1: Section 1. 

Study session 2: Section 2. 

Study session 3: Section 3. 

Study session 4: Section 4 and Section 5. 


Each session requires two to three hours. 

Sections 2 and 3 may take longer to study than the 
other sections. Section 5 is short and will not be 
assessed. 


The optional Video Band C(iv), Algebra 
workout — Differential equations, could be 
viewed at any stage during your study of this 
chapter. 


Introduction 


In many modelling situations, the ultimate mathematical aim is to obtain 
an explicit relationship between two variables. This involves expressing the 
dependent variable, which we denote here by the symbol y, as a function F 
of an independent variable, denoted by x. However, it may not be possible 
to write down an equation of the form y = F() straight away. This desired 
outcome may be achievable only by solving some other type of equation. 


Often the equation that we have to solve involves not only the variables x 
and y themselves but also derivatives of y with respect to x. You saw 
examples of such situations in Chapter C2, Section 3, in the context, of 
modelling motion. Equations that include derivatives, such as 


dy _ Ca ee 
= 2y and aw +6277 + 9x*y = sin(2z), 


are called differential equations. Any function y = F(x) which satisfies 
such an equation is called a solution of the differential equation. For 
example, the function y = e* is a solution of the first differential equation 
displayed above. 


The order of a differential equation is the order of the highest derivative 
that appears in it. So the first equation displayed above is an example of a 
first-order differential equation, since it contains only the first derivative of 
the dependent variable with respect to the independent variable. A 
second-order differential equation, such as the second equation displayed 
aboye, contains a second derivative, and maybe a first derivative as well, 
but no higher derivatives. 


This chapter concentrates exclusively on first-order differential equations, 
and only on particular types of these. Section 1 deals with differential 
equations that can be written in the form 


dy _ 
= f(x), 


where f is a given function. Such equations can be solved by direct 
integration of the function f. Section 2 shows how differential equations 
that can be written in the form 


oY = s(o)aly) 


may be solved, using the method of separation of variables. This method 
again depends on being able to perform integrations. 


Section 3 looks more closely at a special type of differential equation 
soluble by separation of variables, namely dy/dx = Ky, where K is a 
constant. This type of equation can be used to model decay and growth 
processes, such as radioactive decay or the growth of a population. 


Section 4 shows how the computer can be used to plot direction fields, 
which give a visual impression of the overall behaviour of solutions of a 
differential equation. This is followed by a demonstration of how solutions 
of differential equations may be obtained both numerically and graphically. 
using Puler’s method. These approaches apply more widely than those 
considered earlier in the chapter. 


Finally, Section 5 indicates briefly some directions in which the theory of 
differential equations can be further extended. 


In these examples the 
independent variable was t, 
denoting time, while the 
dependent variable was 

v (velocity) or s (position). 
In this chapter we normally 
use Leibniz notation for 
derivatives. This is common 
in the context of differential 
equations. 


Higher-order differential 
equations are defined 
similarly. 


1 Differential equations and direct integration 


Solutions expressed in the 
form of ‘a more general 
equation relating « and y’ are 
discussed in Section 2. 


See Chapter C2, Section 3. 


These rules were introduced 
in Chapter C1, Section 4. 


In Subsection 1.1 we look at first-order differential equations and their 
solutions. In Subsection 1.2 the method of direct integration is used to 
solve a simple type of differential equation. 


1.1 Differential equations and solutions 


We start by drawing together some of the terminology referred to already 
in the Introduction. 


Definitions 

© A differential equation is an equation that relates an 
independent variable, x say, a dependent variable, y say, and one 
or more derivatives of y with respect to x. 

© The order of a differential equation is the order of the highest 
derivative that appears in the equation. A first-order 
differential equation involves the first derivative, dy/dx, and no 
higher derivatives. 

© A solution of a differential equation is a function y = F(x) (or a 
more general equation relating 2 and y) for which the differential 
equation is satisfied. 


As you saw in Chapter C2, a differential equation may feature symbols for 
the independent and dependent variables other than x and y, and you will 
see further examples of this later in the current chapter, when a modelling 
context makes it appropriate. While considering the general theory of 
differential equations, however, we shall stick to the use of x and y. 


A first-order differential equation features the first derivative, dy/dx. The 
differential equation may also involve separate occurrences of x or y, or 
both. Thus each of the following is a first-order differential equation: 
a ay _ 35? 4 9, By = ay. esa . 
dx dx dx dr 3y? +1 
You will shortly see methods which lead to the solutions of all these 
differential equations. First, however, you are asked to verify in a number 
of cases that a function which is claimed to be the solution of a differential 
equation does in fact satisfy that equation. This type of exercise is useful 
in checking that no error has been made in the working leading to a 
ice of checking whether a calculated 
indefinite integral is correct by differentiating the answer obtained. 


In addition, such checking will provide you with some further practice in 
differentiation, using Leibniz notation, including application of the 
Product, Quotient and Composite Rules. The following example indicates 
what is required. 
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Example 1.1 Verifying that a given function is a solution 


Tn each case, show that the given function satisfies the given differential 
equation. 


(a) y=2° +22; dy _ag240 


dx 
2 dy 1-2? 
(b) US Te dx (142) 
dy (3+2xr)y 
= ieee. ka 
(c) y=2*e”; ie ss (x > 0) 


(@) y=inatz4); Y= 220 


Solution 

Tn each case, we differentiate the given function to find an expression 

equivalent to the left-hand side of the differential equation, and then show 

that the same expression arises from the right-hand side of this equation. 

This verifies that the function is a solution of the differential equation. 

(a) The differential equation is dy/da = 327 + 2. The derivative of the 
given function, y = x* + 22, is 


d 3 
=a 2r) = 327 +2, 
a + 22) = 32° + 
which is the same expression as the right-hand side of the differential 


equation, as required, 


The differential equation is dy/da = (1 — x”)/(1 + .2*)*. In order to 
differentiate the given function, y = «/(1 + 2*), we apply the Quotient 
Rule. The function can be written as y = u/v, where 


fe 


u=a and v=1+a2*, 


for which 
du dv 
ee — =22. 
dx andl a 
Hence we have 
dy (ef ez me) 
dr 
_ (1+2?)(1) = ()2r) 
(14+ 27)? 
_l+2?=22? 1-27 
~ (1422)? (1+2?)?' 
which is the same expression as the right-hand side of the differential 
equation, as required. 
(c) The differential equation is dy/dx = (3 + 2x)y/x. In order to 
differentiate the given function, y = r°e?", we apply the Product Rule. 
The function can be written as y = uv, where 


u=2" and v=e™, 


for which 
du_ > Ce 
ae 32° and a 2e’ 


We describe each function F 
by writing y = F(x). 


See Chapter C1, 
Subsection 4.2. 


See Chapter C1, 
Subsection 4.1. 


The restriction « > 0 in the 
statement of part (c) ensures 
that division by « is valid 
here. Another suitable 
restriction is x < 0. 


See Chapter C1, 
Subsection 4.3, 


Recall that expa is an 
alternative (and in this case 
convenient) way of writing e*, 
and that exp is the inverse 
function of In. 


Hence we have 
dy du dv 
ae = rie F LEE 
= (32?)(e2*) + (2*)(2e*) = (3 + 2x) xe”. 
Now we substitute y = r*e** into the right-hand side of the differential 
equation, to obtain 


Sede 
(3 + 2ax)y es (3 + 2r)2*e =(3-+2x)z2e*. 


os x 
This is the same expression as was obtained above for the derivative, 
as required. 


=z 


The differential equation is dy/dx = 2xe~". In order to differentiate 
the given function, y = In(1 + ?), we apply the Composite (Chain) 
Rule. The function can be written as 


y=Inu, where u=1+2", 


for which 

dy 1 du 

duu de 
Hence we have 

dy _ dy du 

dx dudx 

Dove 5 ae. 
=) "ip 


Now we substitute y = In(1 + 2”) into the right-hand side of the 
differential equation, to obtain 
Qae-Y = 2x exp(—In(1 + 2*)) 
= 2rexp(In((1 + 27)7')) 

22 
1+2?° 
This is the same expression as was obtained above for the derivative, 
as required, 


= 22(1 +27) = 


Here are some similar cases for you to try. 


Activity 1.1 Verifying that a given function is a solution 


In each case, show that the given function satisfies the given differential 
equation. 


(a) y=5yz; ONES (x > 0) 
(b) y =sin(e**); 


(c) y = e* cos(2z); u = y(1—2tan(2x)) (-i7 <x < jn) 


i ee ee 2 
@ y= B=-w+y 


Solutions are given on page 43. 
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In Example 1.1(a), we showed that the function y = x° + 2x is a solution 
of the differential equation dy/dx = 3x? + 2 because 


di inet 
a + 2x) = 32? + 2. 
However, this is not the only solution of the differential equation. Since 
d 3 =! 2 
a? +22 +c) = 32" + 2, 


where ¢ is any constant, any function of the form y = x* + 2x +c is also a 
solution of the differential equation dy/dx = 3x? + 2. We say that 

y=«* + 2x +c, where c is an arbitrary constant, is the general solution of 
this differential equation, since it describes the whole infinite family of 
possible solutions. By contrast, the function obtained when a particular 
value is chosen for c, as with the original function y = x* + 2x (where 

e¢ = 0), is called a particular solution of the differential equation. 


You may recall that an arbitrary constant also arises when finding the 
indefinite integral of a function and, as you will see shortly (or may 
already have noticed), there is a very close link between integration and 
the process of solving differential equations. 


The general solution of a first-order differential equation will usually 
feature an arbitrary constant. This does not always occur in the form of a 
+c’ term tagged on at the end of the solution formula, but that will be 
the pattern for the differential equations solved in this section. 


General and particular solutions 


© The general solution of a differential equation is the set of all 
possible solutions of the equation. It usually involves one or more 
arbitrary constants. 


© A particular solution of a differential equation is a single 
solution of the equation, which consists of a relationship between 
the dependent and independent variables that contains no 
arbitrary constant. 


In many cases, where a differential equation arises in a mathematical 
model, the value of the dependent variable is known for one value of the 
independent variable. For example, in Chapter C2 the motion of a 
parachutist was modelled by a differential equation for dv/dt, where v is 
velocity and ¢ is time. In addition, a value of v was specified at time t = 0. 
This permitted an appropriate value to be chosen for the arbitrary 
constant that arose from integration, and hence led to a single function (a 
particular solution) which described the velocity of the parachutist. 


Similar considerations apply more generally. A first-order differential 
equation has a general solution and, in order for a particular solution to be 
obtained from this, a further condition on the solution is needed. This 
often takes the form of an initial condition, in which a value of the 
dependent variable y is specified at a given value of x. When taken 
together, the differential equation and initial condition are called an 
initial-value problem. 


For a discussion of integrals 
and arbitrary constants, see 
Chapter C2, Section 1. 


See Chapter C2, Example 3.1. 
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Initial-value problem 


© An initial condition associated with a first-order differential 
equation requires that the dependent variable y takes a specified 
value, b say, when the independent variable x has a given value, a 
say. This is often written as 


The second notation stems y=bwhenr=a, oras y(a)=b. 
ae eomunion practice of The numbers a and b are called initial values for x and y, 
respectively. 
y= ula), © The combination of a first-order differential equation and an 
in which the symbol y does initial condition is called an initial-value problem. The 
‘double duty’, as a variable on solution of an initial-value problem is a particular solution of the 


the left-hand side and as a 


‘ : bee differential equation which also satisfies the initial condition. 
function on the right. This is 


convenient because it reduces 
the number of separate The following example and activity concern finding a particular solution 


symbols that might otherwise — from the general solution of a differential equation, using an initial 
be required. Nor does it cause eondition 


any ambiguity, since the 
context determines which of 
the two uses is intended. 


Example 1.2 Applying an initial condition 


This is revision of what you As you have seen, the general solution of the differential equation 
saw in Chapter C2. dy/dx = 3x? + 2 is y = 2° + 2x +c, where c is an arbitrary constant. Find 


Example 3.1, Activity 3.1 and the particular solution that satisfies the initial condition 
Subsectior 


y=9whenz=2; thatis, y(2)=9. 
Solution 


We substitute the values 2 = 2 and y = 9 (simultaneously) into the 
expression y = x° + 22 +c. This gives 


9=B42x2+c=12+e, 
so ¢ = —3. The required particular solution is therefore 


y=2°+2r-3. 


Here is a similar problem for you to try. 


Activity 1.2 Applying an initial condition 


The general solution of the differential equation 
dy’ 5 


= 

dr 2/x ( ) 
This can be seen by adapting is y = 5\/x + ¢, where ¢ is an arbitrary constant. Find the particular 
the solution to solution that satisfies the initial condition 


Activity 1.1(a). 
eee) y=Twhenr=1; thatis, y(1)=7. 


A solution is given on page 43. 
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1.2 Direct integration 


In the remainder of this section we concentrate on differential equations 
which can be written in the form 
dy 
ae Tf): (1.1) 
As you saw in Chapter C2, differentiation is ‘undone’ by integration, and 
so the general solution of the differential equation (1.1) is given by the 
indefinite integral of the function f(); that is, 


v= [ f(var. 


Example 1.3 Finding the general and a particular solution” 


(a) Find the general solution of the differential equation 
B32 
dx : 


(b) Find the particular solution of this differential equation that satisfies 
the initial condition y = 2 when x = 1. 


Solution 


(a) The general solution of the differential equation is given by 


y= [ eae. 


The general solution is therefore 
3 


where c is an arbitrary constant. 


(It is always possible, and advisable, to check a general solution by 
differentiating and substituting back into the differential equation. 
Here y = 42° + ¢, so dy/dx = 2?, as required.) 
(b) Putting 2 = 1 and y = 2 into the general solution, we obtain 
2=3(1) +e. 
Hence ¢ = 5 and the required particular solution is 


= 1,545 
yssr +3: 


The method used in Example 1.3 is called the direct integration method 
for solving differential equations. It can be applied to any differential 
equation, such as dy/da = x*, for which the derivative of the dependent 
variable is equal to a known function of the independent variable. 


See Chapter C2, Section 1. 
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Recall that F(a) is an 
integral of f(a) if 

d 

(F(a) = f(a). 
Hence y = F(x) is a 
particular solution of the 
differential equation. 


We say ‘in principle’ because 
in practice it may not be 
possible to find explicitly an 
integral of f. However, you 
will not meet such cases in 
this text. 


You are not asked to solve 
these differential equations. 
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Direct integration 
© The general solution of the differential equation 
dy _ 
Ts f(x) 
is the indefinite integral 
y= [saa =F(z)+¢, 


where F(x) is any integral of f(x) and c is an arbitrary constant. 
© Any initial condition 
y(a) = b, 


enables a value for the arbitrary constant c to be found. The 
corresponding particular solution satisfies both the differential 
equation and the initial condition. 


y=bwhen z=a, that is, 


The method can also be applied, of course, if symbols other than x and y 
are used for the independent and dependent variables. Any differential 
equation of the form 
d(dependent variable) __ ,,. 2 
ieemat 
can in principle be solved using direct integration. For example, the 


differential equation 
dq 
=— =4t+3 
at +3 


has the general solution 
qe [se+syae = 2 43t+e, 


where c is an arbitrary constant. 


Activity 1.3 Can direct integration be applied? 


To which of the following differential equations can the direct integration 
method be applied? 


dy dx _ dz _ 
(a) — =xcosx (b) eT TOT (c) Fe = boost 


dx 


dy ie dx _ da 
(d) ae 7 Noes (e) ge oY (f) ae Tost 


Solutions are given on page 43. 
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Activity 1.4 Applying direct integration 
(a) Find the general solution of the differential equation 
dy _ oa 
ae 
(b) Find the particular solution of this differential equation that satisfies 
the initial condition y = 2 when xz = 0. 
Solutions are given on page 44. 


In Activity 1.4 you saw that the general solution of the differential equation 
dy/dx = e** is y = }e?* +c, where c is an arbitrary constant. Some. graphs 
corresponding to this family of solutions are shown in Figure 1.1. Each of 
these solution curves has the property that, at any point (7, y) on it, the 
gradient of the curve is dy/da = e?". All the solution curves in Figure 1.1 
are parallel, in the sense that the vertical separation of any two curves is a 
constant, independent of x. This is a property of all differential equations 
of the type dy/dx = f(x). However, as you will see in the rest of the 
chapter, it is not a general property of first-order differential equations. 


Activity 1.5 Applying direct integration again 
Solve the initial-value problem 
de... 
ae =sint, 2=1whent=0. 


A solution is given on page 44. 


Activity 1.6 Mass of a burning rocket 


A toy rocket consists of a casing and fuel. Initially the total mass of the 
rocket is 200 grams. The fuel is ignited at time ¢ = 0 and then burns in 
such a way that the rocket’s total mass m grams at time t seconds satisfies 
the differential equation 


(t > 0). 

When the fuel is exhausted, the residual mass of the rocket casing is 
100 grams. 

(a) Find the general solution of the given differential equation. 


(b) Find the particular solution of the differential equation that satisfies 
the initial condition m = 200 when ¢ = 0. 


(c) Hence find the time at which the fuel is exhausted: that is, find the 
time t when m = 100, for the particular solution obtained in part (b). 


Solutions are given on page 44. 


You can check the general 
solution by differentiating it 
and then substituting back 
into the differential equation. 


Figure 1.1 Graphs of 
y= $e* +e, fore =—$,0 


3 
and 3 


In this activity the 
independent variable is ¢ and 
the dependent variable is x. 


The derivative dm/dt 
represents the rate of change 
of mass. This quantity is 
negative because the mass 
decreases as the fuel is burnt. 
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You are asked to establish 
this differential equation in 
Exercise 2.4. 
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Summary of Section 1 


This section has introduced: 

© first-order differential equations, and the solutions of such equations 
(general and particular); 

© the process of obtaining a particular solution from a general solution, 
using an initial condition; 

© the method of direct integration, to solve differential equations of the 
form dy/dx = f(x). 


Exercises for Section 1 


Exercise 1.1 
Show that the function y = —1/zx (x > 0) is a solution of the differential 
equation dy/dx = y?. 


Exercise 1.2 
(a) Find the general solution of the differential equation 
ed =Vr (x>0). 


(b) Find the particular solution of this differential equation that satisfies 
the initial condition y = 5 when x = 4. 


(c) Hence find the value of y when x = 25. 


Exercise 1.3 


Solve each of the following initial-value problems. 


(a) tis cos(2x), u=—2 when x = }7 


dx 

(toe | 
(b) =z (t>0), e=3whent=1 
Exercise 1.4 


A large spherical snowball melts in such a way that, at any instant, the 
rate of decrease of its volume V is proportional to its surface area A, with 
constant of proportionality k > 0. Initially the snowball has a radius of 
100cm. After 2 days the snowball has lost half its initial volume. 


Using the formulas V = $7r* and A = 4rr? for the volume and surface 
area of a sphere, in terms of its radius r, it can be shown that the change 
in radius r (in cm) of the snowball with time ¢ (in days since melting 
started) is modelled by the differential equation 


(a) Find the general solution of this differential equation. 


(b) Find the particular solution of the differential equation that satisfies 
the given initial condition. 


(c) Use the condition that the snowball has lost half its initial volume 
after 2 days to find the value of the constant of proportionality k. 


(d) How long will it take for the snowball to disappear completely? 


2 Solution by separation of variables 


In Subsection 1.2 you saw that the general solution of the differential 
equation 

dy 

ee f(x) 


is given by the indefinite integral 
y= [ f(vax. 


In this section we develop a method for finding solutions of a further type 
of differential equation, namely, those of the form 


dy 


de = 19H), 


where f and g are known functions. This method, called separation of 
variables, is described in Subsection 2.2. As a preliminary, Subsection 2.1 
shows how the Composite (Chain) Rule for differentiation may usefully be 
employed even when one of the component functions in the composite is 
not known explicitly. 


2.1 Further use for the Chain Rule 


In Section 1 you saw the solutions of several differential equations. In each 
case, the solution was specified in the form y = F(z), for some function F. 
When solving a differential equation. it is sometimes not possible to reach 
such an explicit form for the solution function. For example, suppose that 
it has been deduced, from a differential equation for dy/dx, that the 
variables x and y are related by the equation 


yty=2?+o. (2.1) 


It is not straightforward to solve this equation for y and hence to reach the 


form y = F(a) for the solution. However, equation (2.1) does describe a The possibility of a solution 
solution, even if only implicitly. For example, a corresponding solution being ‘a more general 
curve can be plotted, and the value of dy/dz at a point (x,y) on this curve equation relating x and y’ was 
is the slope of the curve at that point. allowed for in the definition of 
a solution of a differential 

More generally, suppose that the solution of a differential equation is equation in the box on page 6. 
obtained in the form 

H(y) = F(2), (2.2) 


where H and F are known functions, and H(y) 4 y. (Equation (2.1) is of 
this form.) This is called an implicit solution of the differential equation. 
By contrast, a solution written in the form y = F(x), where F is a known 
function, is called the explicit solution of the differential equation. In 
seeking to solve a differential equation, it is the explicit solution which we 
seek. However. if this is impossible to find, then we must settle for an 
implicit solution. 


In Subsection 1.1 you saw how to verify that a given function is an explicit See Example 1.1 and 
solution of a specified differential equation. In order to check that a Activity 1.1. 
calculated implicit solution satisfies the differential equation from which it 

was obtained, the Chain Rule is employed. This is possible because, in 

equation (2.2), the y that appears on the left-hand side is (implicitly) a 


15 


CHAPTER C3 DIFFERENTIAL EQUATIONS AND MODELLING 


function of x. (All the implicit solutions we deal with have the form of 
equation (2.2) and are such that y is a function of x.) If we write 


Here = is an additional z= H(y) = F(z), 
variable, introduced only to 


assist the ‘argument at thia then two expressions can be deduced for dz/dz as follows. Since z = F(x), 


point. we have 
dz 
SSF 
ie (x) 
Since also z = H(y), where y is a function of x, we have, by the Chain 
Rule, 
dz _ de dy _ py dv 
dx dydx dx” 
(This step is called implicit differentiation of H(y) with respect to x.) 
Equating these two expressions for dz/d:r, we see that 
di 
if H(y) = F(x), then Hy) = F"(z). (2.3) 
The next example indicates how this result can be put to use. 
Example 2.1 Verifying a solution given implicitly 
(a) Show that equation (2.1). 
yity=r°+9, 
is an implicit solution of the differential equation 
dy __ 2@ 
dx 3y?+1° 
Note that (1,2) does lie on (b) What is the slope of the corresponding solution curve at the 
the given curve, since point (1,2)? 
24+2=10=17+9. Solution 


(a) We have H(y) = F(x), where H(y) = y+ y and F(x) = 2? +9. It 
follows from result (2.3) that 


which leads to 
(3y? + yw = 22. 
Dividing through by 3y* + 1 gives 
» dy 2a 
de 3 +1" 
as required. 
(b) At (1,2), the slope of the curve is 
dy 2x1 2 


dr 3x41 13° 


Here are similar problems for you to try. 
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SECTION 2 SOLUTION BY SEPARATION OF VARIABLES 


Activity 2.1 Verifying a solution given implicitly 


(a) Show that the equation 
ytsny=2+e*—1 
is an implicit solution of the differential equation 
ote (—17<y<7)- 
(b) What is the slope of the corresponding solution curve at the 
point (0,0)? 
Solutions are given on page 44. 


Activity 2.2 Another verification 
In Chapter Cl you saw that 
d 2 
qa tan) = sec 2. 


Use this result, and the trigonometric identity sec? @ = 1 + tan*@, to show 
that the equation 


tany =r 
is an implicit solution of the differential equation 
dy 1 1 1 
=—=——— (-39<y< 42). 
ie iter Oat <2 <an) 


A solution is given on page 44. 


Comment 


In fact, it is possible in this case to write the solution explicitly. For if 
tany =a and —}7<y< 3m, then y = arctan. Hence the result of this 
activity shows that 


d 
ay erctanz) = 


1+27° 
It can be shown in a similar way that 
1 


Le essing (-1<2<]1). 
dx 


The previous example and activities indicate the usefulness of implicit 
differentiation using the Chain Rule, as described by result (2.3). In the 
next subsection, you will see that result (2.3) enables us to explain a 
solution process for a wider class of differential equations than those 
considered in Subsection 1.2. 


See Chapter C1, 

Activity 4.4(a). The same 
differentiation appears in 
Chapter C2, Activity 2.4(a). 


This trigonometric identity 
was introduced in 
Chapter C2, Subsection 2.1. 


This derivative formula is 
valid for all values of x. 
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The phrase ‘in implicit form’ 
is used to describe a general 
or particular solution that is 
an implicit solution, The 
phrase ‘in explicit form’ is 
applied similarly. 


This demonstrates that. when 
both sides of an equation are 
integrated, we need introduce 
only one constant of 
integration. We shall follow 
this practice from here on. 
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CHAPTER C3 DIFFERENTIAL EQUATIONS AND MODELLING 


2.2 Separation of variables 


We now derive a method for solving differential equations of the form 


where f and g are known functions. The starting point for explaining this 
method is the result, obtained in Subsection 2.1, that 


Hy)! 


y 7 
— = F(x). 
a (x) 
Integration is the reverse process to differentiation. Hence, faced with a 
differential equation of the form 


d 
hy) = f(z). 


where f and h are known functions, a solution will be given by the 
equation H(y) = F(x), provided that 


H'(y)=h(y) and F(x) = f(x). 
It follows that 
d 
if hy) = f(x), then if Aly) ey = | f(x) de. 


The latter equation gives the general solution of the differential equation, 
in implicit form. 


if H(y)= F(x), then (2.3) 


(2.4) 


The next example shows how this result may be applied. 


Example 2.2 Finding a general solution in implicit form 
(a) Find, in implicit form, the general solution of the differential equation 


dy _ 2 
25, F- 


(b) Find the corresponding particular solution that satisfies the initial 
condition y = 3 when x = 2. 
Solution 


(a) Comparing the differential equation with result (2.4), we have 
h(y) = 2y and f(x) = 3x?. Hence the general solution of the 
differential equation is given by 


[vay = [32 de; that is, y2+a=2°+5, 
where a and b are arbitrary constants. These two constants can be 


combined to form the single arbitrary constant ¢ = b— a. Thus the 
general solution, in implicit form, is 


y=2 +e, 
where c is an arbitrary constant. 
(b 


To apply the initial condition, put x = 2 and y = 3 into the general 
solution. This gives 


P=PB+e, 
so c= 1. Hence the required particular solution is 


y=rtl 


SECTION 2 SOLUTION BY SEPARATION OF VARIABLES 


Here is a similar problem for you to try. 


Activity 2.3 Finding a general solution 


(a) Find, in implicit form, the general solution of the differential equation 
adv 
y? dx 

(b) Find the corresponding explicit form of this general solution. 


-1 (y>0). 


Solutions are given on page 45. 


The approach demonstrated in Example 2.2 and Activity 2.3 can be 
extended to solve any differential equation of the form 


Gs ap. re 
dg 1 @)at)- (2.5) 
A differential equation of this type is said to be separable, and the 
corresponding method of solution is called separation of variables. To 
solve equation (2.5), we first rewrite it so that there is a function of x alone 
on the right-hand side. This is done by dividing both sides by g(y), to give 
1 dy 
gly) de dla) 


Now we can apply result (2.4), with h(y) = 1/g(y). This shows that the 
general solution of equation (2.5), in implicit form, is 


[aqua [sorae. (2.6) 


To summarise, the method of separation of variables is as follows. 


Separation of variables 
The method applies to differential equations of the form 


dy . 
=, = F()9(y)- (2.5) 
1. Divide both sides by g(y), to obtain We assume here that any 
jen value y for which g(y) = 0 is 
pee = hth excluded from consideration. 
gly) dx 
2. Integrate both sides with respect to x. From result (2.4), the 
outcome is 
[a dy = [1@ dr. (2.6) At this stage the variables are 
gu) separated, with a function of 
3. Carry out the two integrations, introducing one arbitrary y alone on the left-hand side 
constant, to obtain the general solution in implicit form. If and a function of « alone on 


possible, manipulate the resulting equation to make y the subject, the right-hand side. 


thus expressing the general solution in explicit form. 
(A general solution, in implicit or explicit form, can be checked by 
substituting back into the differential equation.) 
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The combination of a 
differential equation and an 
initial condition is called an 
initial-value problem — see 
page 10. 


Some differential equations of 
the type given by 

equation (2.7) will be 
considered in Section 3. 


CHAPTER C3 DIFFERENTIAL EQUATIONS AND MODELLING 


Any initial condition provides a value for the arbitrary constant in step 3. 
The corresponding particular solution satisfies both the differential 
equation and the initial condition. 


You might like to note that the equation in step 2 can be written down 
quickly from the differential equation obtained in step 1 by regarding 
dy/dx temporarily as an ordinary quotient and ‘multiplying through by dx’ 
before adding the two integral signs. This shortcut is a convenient aid in 
recalling this step of the method. However, remember that the derivative 
dy/dx is not dy divided by dx. We have not assigned a meaning to either 
dx or dy on their own; they make sense only when combined in a derivative 
such as dy/dx, or when preceded by an integral sign as in [... da. 

This method can be applied also when independent and dependent 
variables other than x and y are used, provided that the derivative is equal 
to a product of a function of the independent variable and a function of 
the dependent variable. 


Activity 2.4 Can separation of variables be applied? 


To which of the following differential equations can the separation of 
variables method be applied; that is, which of them are of the form 

d(dependent variable) 2 p , 
Asie ee ble I 
GGadepandentacabial f (independent variable)g(dependent variable) 
for some choice of the functions f and g? For each of the equations below 
that is of this form, identify the functions f and g. 


dy _ 4 aes ») oP 
(a) Fiehied | (b) it =tcost (ec) qe = PCOS 
do mae eS c/a 
(d) oF =ut+y (e) a u+vy (f) de = cos(ry) 


Solutions are given on page 45. 


Parts (b) and (c) of Activity 2.4 illustrate that both the equations 


dy _ 

Pie f(x) (1.1) 
and 

d. 

= = gly) ex 


are special cases of the equation dy/dx = f(x)g(y), with g(y) = 1 in the 
first case and f(x) = 1 in the second. As you saw in Subsection 1.2. the 
general solution of equation (1.1) can be found by direct integration, to 
which separation of variables reduces when g(y) = 1. 


Example 2.3 Applying separation of variables 
(a) Find the general solution of the differential equation 
a =e “cosx. 


(b) Find the particular solution of this differential equation that satisfies 
the initial condition y = 1 when x = 0. 


SECTION 2 SOLUTION BY SEPARATION OF VARIABLES 


Solution 
(a) The differential equation is of the form of equation (2.5) with 


(b 


f(x) =cosx and g(y) =e’. Dividing both sides of the differential 
equation by e~" gives 


d, 
oF = cos x. 


From this we obtain, using result (2.4), 


[eray= [cosede. 


On performing the integrations and including a single arbitrary 
constant c, this becomes 


e’ =sinr +c, 
which is an implicit form of the general solution. 


To make y the subject of this equation, recall that In(e”) = y. Hence, 
on taking the natural logarithm of each side of the equation, we obtain 


y = In(sinz +e), 
where c is an arbitrary constant. This is the explicit form of the 
general solution. 
Putting « = 0 and y = 1 into the explicit form of the general solution, 
we have 

1 = In(sin0 + ¢) = Ine, 
e, and the required particular solution is 


Hence ¢ = 


y = In(e +sinz). 


Comment 


a 


2. 


The solution could now be checked, though the details are not given 
here. 

In part (b) the explicit form of the general solution was used in finding 
the required particular solution. This particular solution can also be 
obtained by using the implicit general solution. 


Putting « = 0 and y = 1 into 
e’ =sine +e, 
we obtain 
e=sin0+e, 
so ¢ =e. Hence the particular solution in implicit form is 


eY=e+sina, 


and taking the natural logarithm of each side gives 
y =In(e+sinz), 
as before. 


Unless you are directed otherwise, you should find particular solutions 
from the explicit form of the general solution. 


Other choices for f(r) and 
g(y) are possible, for example, 


f(z) = 


choices lead to the same 
general solution. 


Since 1/eY =e", dividing by 
e~” is the same as 
multiplying by e”. 


Recall that In(e) = 1. 
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Recall that sec x = 1/(cos a). 
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CHAPTER C3 DIFFERENTIAL EQUATIONS AND MODELLING 


Here are some similar problems for you to try. 


Activity 2.5 Applying separation of variables 
(a) Find the general solution of the differential equation 
dy 
dx 
(b) Find the particular solution of this differential equation that satisfies 
the initial condition y = 1 when x = 0. 


| w>0). 


Solutions are given on page 45. 


Activity 2.6 Applying separation of variables again 
(a) Find the general solution of the differential equation 


dx 

ae 7 secre (-30 <x < }n). 

(b) Find the particular solution of this differential equation that satisfies 
the initial condition 2 = 0 when t = 0. 


Solutions are given on page 46. 


A possible snag in applying the separation of variables method was 
mentioned in a margin note on page 19. We can divide both sides of the 
equation 


oY = F(2)g() 
ar 


by g(y) only if g(y) 4 0, since division by zero is an undefined operation. 
It is therefore necessary to exclude from consideration any value of y, the 
dependent variable, for which g(y) = 0. This may be done by imposing a 
suitable condition on y from the outset, as we did in Activity 2.5 by 
assuming that y > 0. 


As well as having to impose restrictions on the dependent variable in order 
to apply the separation of variables method, we must be aware of other 
pitfalls that can arise. The solution obtained may be valid only for a 
certain range of values of the independent variable. For example, in 
Activity 2.5 you found that the required particular solution of the given 
differential equation dy/dx = —2x/y, assuming that y > 0, was 


y=Vv1—-2. 


The expression y1 — x? is a real number only for values of « such that 

-—1 <2 <1. Since x = 1 and x = =1 both give y = 0, which is excluded by 
the restriction y > 0 in the question, the domain of this particular solution 
is the interval (—1,1). Similar considerations apply to the corresponding 
general solution in Activity 2.5, which was 


y=vb—2?. 


This formula provides positive values for y only when the arbitrary 
constant b is positive, and then only when —b'/? < x < b!/*. So the 
interval (—b!/?,b!/?), where b > 0, is the domain of this solution. 


SECTION 2 SOLUTION BY SEPARATION OF VARIABLES 


These examples show that consideration of the appropriate domain for a 
solution can become a bit complicated. This should not be a cause for 
your concern, as you will not be expected to derive the appropriate domain 
for a solution in the manner indicated above. However, it is as well to be 
aware that in practice the choice of domain may require some thought. 
Deriving a formula y = F(x) (or H(y) = F(:)) for a solution does not 
guarantee that the formula obtained makes sense for all real numbers <r, 
nor that the function F actually satisfies the differential equation without 
some restriction being imposed on its domain, Also, as has been pointed 
out several times, it is wise to check a proposed solution by substituting it 
into the differential equation. 

The final activity of this section investigates a model for the motion of a 
space probe launched from Earth. The modelling theme is continued in the 
next section, where each model features a differential equation that can be 
solved by separation of variables. 


Activity 2.7 Motion of a space probe 


A motorless space probe is launched vertically upwards from the Earth's 
surface. When the probe is at a distance r from the centre of the Earth, its 
outward velocity v satisfies the differential equation 
dy sek 
dr ur? 
where R is the radius of the Earth and k is a positive constant. All the 
variables and constants are measured in the appropriate SI units. 
(a) Motion of the probe outwards from the Earth corresponds to v > 0. 
Find the general solution of the differential equation in this case. 
(b) The initial condition corresponding to a launch velocity of up is 


(r> R), 


v(R) = vp: that is, v =vp_ when r= R. 


Show that the corresponding particular solution for the motion 
outwards from Earth is 


v= /uR +2 ( 


Show that at very large distances from the Earth, the velocity of the 
space probe tends towards the limiting value 


{ec 


u-=- 


R 


Hence, by considering the cases vg > \/2k/R and ug < \/2k/R, 
explain the significance of the particular launch velocity 


2k 
UR=IR- 


Solutions are given on page 46. 


The restriction r > FR arises 
because the differential 
equation is a suitable model 
only when the space probe is 
above the surface of the 
Earth. 


Do not spend too long on 
part (c) if you find it difficult. 
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Summary of Section 2 


This section has introduced: 


© the contrast between the explicit solution of a differential equation, in 
the form y = F(x), and an implicit solution, in the form H(y) = F(x); 


© implicit differentiation of a function H(y) with respect to x, where y is 
a function of x whose rule is not known explicitly; 


© the method of separation of variables, to solve differential equations of 
the form dy/dx = f(x)g(y)- 


Exercises for Section 2 


Exercise 2.1 
(a) Show that the equation 
y+Iny = (1 —2z)e™ 
is an implicit solution of the differential equation 
dy _ y(1—6r)e* 
dx yt+l yd): 


(b) What is the slope of the corresponding solution curve at the 
point (0.1)? 

Exercise 2.2 

(a) Find the general solution of the differential equation 


od 
uw? = eas 
(b) Find the particular solution of this differential equation that satisfies 


the initial condition u = 1 when x = 1. 


Exercise 2.3 


Solve in explicit form each of the following initial-value problems. 


dy _ os eae oi 
(a) Fela (y>0), y= 3 whenr=0 
(b) ate (c>0), 2=2whent=0 
Exercise 2.4 
In this exercise you are asked to establish the differential equation 
dr 
S=-k (r>0 
a 20) 


which was solved in Exercise 1.4. The equation models the change in 
radius r (in em) of a snowball with time ¢ (in days since melting started). 


The snowball melts in such a way that, at any instant, the rate of decrease 
of its volume V = 4zr* is proportional to its surface area A = 4zr?, with 
constant of proportionality k > 0. Use these formulas for V and A, and the 
Chain Rule, to establish the differential equation. 


SECTION 2 SOLUTION BY SEPARATION OF VARIABLES 


Exercise 2.5 


This exercise concerns a mathematical model for the flow of water from a One application of this model 
tank, of constant horizontal cross-section, through a small hole near its is to the water emptying from 
base. According to the model, the height h of the water surface above the —_@ vertically-sided bath via the 
hole is related to the time ¢ since the hole was uncovered. This relationship Plug-hole. 
is described by the differential equation 
dh _ 
dt 
where k is a positive constant. (The value of k depends on both the area of 
the hole and the horizontal cross-sectional area of the water column in the 
tank.) 
(a) Find the general solution of this differential equation. 
(b) Find the particular solution for which h = hp at t = 0. (Here ho is the 
height of the water column above the hole at the moment when the 
hole is uncovered.) 


kh? (h>0), 


(c) Show that, according to this model, the tank empties to the level of 


the hole after a time 2h}/"/k. 
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3 Modelling growth and decay 


In Chapter B3, 

Subsection 4,1, a particle was 
defined as a material object 
whose size and internal 
structure may be neglected, 
The usage here is taken from 
physics, where ‘particle’ 
means an object which is 
finite but very small, 


If time ¢ is measured in days, 
then the unit of k is day~!. 
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In this section we investigate mathematical models that are based on 
first-order differential equations of the form 

dy 

— = Ky, 

dx: _ 
where K is a constant. This type of equation is soluble by separation of 
variables. The resulting behaviour is that of exponential growth or decay. 


Subsection 3.1 concerns radioactive decay, while Subsection 3.2 returns to 
the topic of population change. 


3.1 Radioactive decay 


As you may know, the nuclei of radioactive substances disintegrate 
spontaneously, at a rate which varies according to the substance concerned. 
For example, a gram of pure uranium is composed of an enormous number 
(about 2.5 x 107") of uranium atoms. Each of these atoms consists of a 
central nucleus surrounded by a cloud of 92 electrons. There are three 
different kinds of naturally-occurring uranium nuclei, called uranium-234, 
uranium-235 and uranium-238, but uranium as mined is more than 99% 
uranium—238. (The numbers 234, 235 and 238 give the atomic mass of 
these three kinds of uranium.) From time to time one of the uranium-238 
nuclei disintegrates, emitting an electrically charged particle, called an 
alpha particle, which can be detected using a Geiger counter. As a result: 
of these disintegrations, the uranium slowly decays into another substance, 
namely, thorium-234. The amount of uranium contained in a lump of 
matter therefore decreases steadily with time. 


Each nucleus of uranium-238 has a known probability of decay in a given 
time interval. However, because of the enormous number of nuclei present 
in even one gram of uranium~238, it is sensible to use a continuous model 
to describe this situation. measuring the amount of uranium in grams, 
rather than using a discrete model, where we would measure the number of 
uranium nuclei. The same applies to the decay of other radioactive 
substances. 


As each nucleus has a constant probability of decay in a given time 
interval. we assume for the continuous model that the rate of decay of a 
radioactive substance is proportional to the amount of that substance 
which is present. Consequently, the change in mass m of the radioactive 
substance that is present at time ¢ can be modelled by the first-order 
differential equation 

dm 

dt 
where k is a positive constant, called the decay constant. The negative 
sign in this differential equation reflects that we have a decaying process, 
where the mass m is a decreasing function of time, so that dm/dt is 
negative. The condition m > 0 corresponds to the fact that decay can 
occur only when some of the radioactive substance is present. Physically, 
we would expect equation (3.1) to be a good model for the situation. 
except where there are only a small number of radioactive nuclei, for which 
a discrete model would be more appropriate. 


=-km (m>O0), (3.1) 


SECTION 3 MODELLING GROWTH AND DECAY 


To find a particular solution of equation (3.1), we require an initial 
condition. If there is an amount mp of the radioactive substance present at 
time t = 0, then the initial condition is 

m=mp whent=0; that is, m(0) = mo. (3.2) 


The next activity asks you to solve the initial-value problem given by 
equations (3.1) and (3.2). 


Activity 3.1 Solving the radioactive decay problem 


(a) Use separation of variables to find the general solution of 
equation (3.1), giving your answer in explicit form. 

(b) Find the corresponding particular solution which satisfies the initial 
condition (3.2). 

Solutions are given on page 46. 


In the above activity you showed that the initial-value problem that 


models radioactive decay, 
di 
> =—-km (m> 0), m = mp when t = 0. 


has solution 


m = moe"™" 


The graph of this solution, shown in Figure 3.1, illustrates that the 
solution has the characteristics we would expect of radioactive decay. The 
mass of the substance decreases, initially quite quickly and then more 
slowly, until eventually there is practically none of the substance left. 


Different radioactive substances have different values for the decay 
constant k. The larger the value of k, the faster the substance decays. 
This is illustrated in Figure 3.2. 


The exponential decrease of a radioactive substance has an interesting 
consequence in terms of the factor by which the amount of the substance 
declines over any fixed time interval. If we denote the mass of the 
radioactive substance at time ¢ by m(t) then, according to equation (3.3), 
we have 

m/(t) = moe". 
After a further time interval of duration T, the mass of the radioactive 
substance has become 


m(t + T) = moe *t*T) = mge~**e“*T = m(t)e“*7. 
Tt follows that 
m(t+T) _ .-er (3.4) 


mi(t) 
This equation says that, at the end of a time interval of duration T, the 


proportion of the original amount of a radioactive substance remaining is 
e~*T, regardless of the time ¢ at the start of the interval. 


According to equation (3.4). the time which it takes for the mass of 
radioactive substance to diminish to half its original amount is given by 


1 m(t+T) _ er 


2 m(t) 


ma 
mo 
m= moe aoe 
ee 
0 t 


Figure 3.1 Graph of 
m=moe", k > 0 


ma, 


0 1 


Figure 3.2 Graph of 
m = moe~** for several 
values of k 


Here we use the rule 


erty = eel, 
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Recall that 
In} =In(27) = —In2. 


You do not need to be able to 
follow the detail of the 
following outline description 
in order to attempt 

Activit; 3. 


Tn terms of numbers of atoms, 
this proportion is 1 to 10!. 


This discovery was made in 
the 1950s by the American 
chemist Willard F. Libby, He 
was awarded the Nobel Prize 
for Chemistry in 1960 for his 
pioneering work in 
radiocarbon dating. 
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Taking the logarithm of each side, we obtain 


In} =-AT; thatis, kT =In2. 


Hence the mass of radioactive substance decays to half of its original 
amount in the time 
In2 

T=, 
where k is the decay constant. This time T is called the half-life of the 
substance. Whereas the speed of decay of a radioactive substance could be 
described in terms of the decay constant k, it is more usual to use the 
half-life 7. Radioactive substances have a very wide range of half-lives, 
from millions of years down to tiny fractions of a second. 


(3.5) 


Activity 3.2 Using the half-life 


Silicon-31 has a half-life of 2.62 hours. Find the value of the decay 
constant k for silicon-31. Hence calculate the proportion of a sample of 
silicon-31 which will still be present after six hours. Give each value to 
three significant figures. 


A solution is given on page 47. 


Carbon dating 


One useful application of radioactivity is to the dating of archaeological 
finds of biological origin. Carbon occurs naturally in three forms: 
carbon-12, carbon-13 and carbon-14. Chemically they are 
indistinguishable, but carbon-14 is radioactive, with a half-life of 
approximately 5570 years, and decays continually (into nitrogen). This 
decay is balanced by the production of carbon~14 in the atmosphere, and 
it may be assumed that the proportion of carbon-14 in naturally-occurring 
carbon is constant over time. This same known proportion exists in all 
living tissue. 


When an organism dies and is buried, the carbon-1d is no longer 
replenished and the proportion of carbon-14 in the carbon decreases with 
time. The quantity of carbon~14 in a sample of dead tissue can be 
estimated from a measurement of the radioactivity, and the total quantity 
of carbon in the sample can be estimated from chemical analysis. These 
values, the known half-life of carbon—14 and the known proportion of 
carbon-14 in naturally-occurring carbon in living tissue form the basis for 
a method of estimating when the organism died. The method is known as 
carbon dating. In particular, an estimate of the proportion of the 
original amount of carbon-14 still present in the sample can be found. 


Activity 3.3 Using carbon dating 


The bones of an animal are found in an archaeological dig. Analysis 
produces the estimate that 85% of the original amount of carbon-14 is still 
present in the bones. Taking the half-life of carbon—14 to be 5570 years, 
find the approximate age of these bones to three significant figures. 


A solution is given on page 47. 


Finding the decay constant from data 


You have seen that the amount m of a radioactive substance present in a 
lump of matter at time t¢ is modelled by an initial-value problem whose 
solution is 


m=moe™. (3.3) 


The decay constant k cannot be measured directly, but it can be deduced 
by measuring experimentally how the amount of the radioactive substance 
changes with time. The data in the table below were obtained by Ernest 
Rutherford early in the twentieth century, when he was studying 
radioactive decay. He observed a sample of a radioactive compound of 
thorium over a period of several months, taking measurements from time 
to time, and deducing the proportions m/mo of the original mass of the 
sample which remained. = 


48 66 89 
0.22 0.13 0.06 


t(days)| 10 15 20 26 33 


m/mp | 0.71 0.59 0.51 043 0.33 


According to the model, the amount of radioactive substance present, as a 
proportion of its initial mass, is given by m/mo = e~*". Taking the 
logarithm of each side of this equation. we obtain 


In (=) = —kt. 
mo 


Hence a plot of In(m/mp) against t, for the data above, should result in a 
straight line through the origin with slope —k. After adding to the table 
above a row for the values of In(m/mp) to two decimal places, we obtain 
the table below. 


t (days) | 10 15 20 26 33 48 66 89 
m/mo 0.71 059 0.51 043 0.33 0.22 0.13 0.06 
In(m/mo) | —0.84 -0.53 -0.67 —0.84 -111 -1.51 —2.04 —-2.81 


In Figure 3.3 a graph of In(m/mp) against t has been plotted. It shows 
both the data points and a straight line through the origin which fits them 
well. The line is drawn through the origin because, when t = 0, we have 
In(m/mg) = In(mo/me) = In1 = 0. 


In(m/mo). 20 40 60 80 


vn 


Figure 3.3 Graph of In(m/mp) against t for the data in the table above 


mg is the original amount of 
the substance. 


Ernest Rutherford 
(1871-1937) was a New 
Zealand-born physicist of 
great distinction. He was the 
first to show that 

s caused by the 
spontaneous decay of heavy 
atoms into slightly lighter 
ones, for which he was 
awarded the Nobel Prize for 
Chemistry in 1908. He also 
propounded the nuclear 
model for the atom, and 
became the first person to 
split the atom. 


A straight line through the 
origin in the (ic, y)-plane, with 
slope —k, has the equation 

uy = —kax. Here we have ¢ in 
place of x and In(m/mo) in 
place of y. 


In Block D you will see an 
objective method for 
calculating the line which 
best fits a set of data points. 
However, it is often sufficient 
to use one’s eye. This is a 
rather subjective approach, 
but a good rule of thumb is to 
draw the line so as to 
minimise the distances of the 
data points from the line, 
while ensuring that there are 
about the same number of 
data points on either side of 
the line. This rule of thumb is 
employed in Figure 3.3. 


From the graph you could 
read the second coordinate as 
slightly different from —3.2. 
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Figure 3.3 shows that these data points are reasonably well-fitted by a 
straight line through the origin with negative slope. This validates use of 
the exponential model to describe radioactive decay. The straight line 
passes through the origin and the point (100, —3.2), so its slope is 
approximately —3.2/100 = —0.032. Hence the decay constant for 
radioactive decay of the thorium compound is estimated to be 


k = 0.032 day~'. 
From equation (3.5), the corresponding half-life is 


_In2_ In2 
~ k 0.032 

A number of remarks about the above estimation, the decay constant and 

the half-life are called for. 

1. It is assumed in the above method for estimating the decay constant 
that hand plotting of the data and the graph are used. 

2. The values of In(m/mo) were given to two decimal places because hand 
plotting on a normal-sized piece of graph paper can represent this level 
of accuracy at most. 

3. The first coordinate of the point chosen in order to calculate the slope 
of the line drawn by eye should be near the high end of the t data 
range, and ideally be arithmetically convenient. The first coordinate of 
the point (100, —3.2) satisfies both requirements. 

4. The value of T was given to two significant figures. the same accuracy 
as for the value of k, estimated from the graph. Slightly different 
estimates for k, obtained from other lines “fitted by eye’. can affect the 
value of T. For example, the point (100, —3.25) leads to k = 0.0325 
and T = 21.3 to three significant figures. 


= 22 days (to 2 s.f.). 


You should bear the above points in mind when attempting the following 
activity, and in similar situations. 


Activity 3.4 Finding the half-life of 


The table below shows data for the radioactive decay of a sample of 
uranium-239. 


t(minutes)| 10 20 30 40 50 
m/mo 0.73 0.56 0.42 0.29 0.24 
In(m/mo) 


(a) Calculate to two decimal places the entries for the third row of this 
table. 

(b) Plot, on graph paper, a graph of In(m/mo) against t for these data. 
Draw ‘by eye’ a line through the origin which you think best fits the 
data points. 

(c) Using your straight line from part (b), estimate values for the decay 
constant k and the half-life T for uranium-239. 


Solutions are given on page 47. 
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Comment 


As was said before, fitting a straight line to data points ‘by eye’ is a 
subjective process. This means that your values of the decay constant and 
half-life may not be exactly the same as those given in the solution. 
However, they should not be very different! 


3.2 Continuous model for population growth 


Earlier in the course you saw mathematical models, based on recurrence 
relations, to describe the way in which population size alters with time. 
These models all had the property that they predicted population size only 
at discrete intervals of time, typically at annual intervals. The aim here is 
to develop a model which describes a population at any instant of time, on 
a continuously changing basis. In practice, of course, populations are made 
up of individuals, and so the size of a population can take only integer 
values and does not change continuously. However, if the population size is 
large, then its behaviour may be approximated satisfactorily by a 
continuous function, and continuous models for populations are therefore 
of value. 


If P represents the size of a population at time ¢, then the instantaneous 
rate of change of the population size is given by the derivative dP/dt. A 
simple modelling assumption is that the rate of change of the population 
size is proportional to the current population size. We then have the 
differential equation 


dP a 

ars KP (P>0), 

where K is a constant, called the proportionate growth rate. The 
condition P > 0 corresponds to the fact that population change can occur 


only for a positive population. 


The value of K may be positive (birth rate exceeds death rate), negative 
(birth rate less than death rate) or zero (birth and death rates equal). If 
the initial population size (at time ¢ = 0) is Py, then we have the 
initial-value problem 


dP ” 

Pee KP 
This problem has the same mathematical form as that given by 
equations (3.1) and (3.2) for the radioactive decay model. Hence the 
solution to problem (3.6) can be written down immediately. by comparison 
with the previous case whose solution is given by equation (3.3). The 
solution to problem (3.6) is therefore 

P= Pek, 
Figure 3.4 shows some typical graphs of this solution. The case K > 0 
corresponds to population growth, K = 0 gives a static population size and 
K < 0 corresponds to population decline. The behaviour for K < 0 is 
mathematically identical to that examined earlier for radioactive decay 
(where A :), so we shall now look at examples for which K > 0. The 
next activity shows how equation (3.7) can be put to use. 


(P>0), P=P)whent=0. (3.6) 


(3.7) 


Discrete models for 
population change were 
developed in Chapter B1, 
Sections 2-4. 


The only difference is that —k 
had to be negative (k > 0), 
whereas KC (in place of —k) 
can take any real value. 


Figure 3.4 Graph of 
P = Pye** for several 
values of K 
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See Chapter B1, 
Subsection 2.2. 


See Activity 2.2 of 
Chapter B1. See also 
Figure 2.6 of that chapter. 
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Activity 3.5 Exponential growth of a human population 


The population size of a country is currently 5 million, and its 
proportionate growth rate is 3% per year, so that K = 0.03 year~'. 


(a) Use equation (3.7) to predict what the population size of the country 
will be after ten years and after a hundred years. 


(b) After how long will the population size reach 50 million? 
In each part, give your answers to three significant figures. 


Solutions are given on page 48. 


For obvious reasons, this model for population change is called the 
exponential model. You may recall that in Chapter B1 this name was 
given to the first discrete model for population change. Which of the two 
models is being used should be clear from the context, but if there is any 
possibility of confusion, then we distinguish between them by referring to 
either the discrete exponential model or the continuous exponential model. 


As already noted, the continuous exponential model for populations has 
the same mathematical form as that of the radioactive decay model. It 
follows that, in place of equation (3.4), we now have 


PST) _ 
Poe 
is equation says that, in a time interval of duration 7’, the population 
e changes by a factor e*?, regardless of the time t at the start of the 
interval. If K <0, then the population is in decline and has a half-life T, 
given by T = —(In2)/K, just as for radioactive decay. If K > 0, then the 
population size increases. Here we can define the doubling time T for the 
population, which is the time that it takes for the population to double in 
size from any starting value. 


(3.8) 


Activity 3.6 Doubling time 


(a) Suppose that the growth of a population is described by a (continuous) 
exponential model, with proportionate growth rate K > 0. Determine, 
in terms of A’, the doubling time for the population. 

(b) Find the doubling time, to three significant figures, for the population 
described in Activity 3.5. 


Solutions are given on page 48. 


Finding the parameters of the exponential model 


As with the decay constant for radioactive decay, the proportionate growth 
rate for a population may be deduced from data on population sizes. In 
the case of a human population, this will be census data. 


For example, the table below shows census data for the USA for the years 
1790-1890. We analysed these data to an extent in Chapter B1, using the 
discrete exponential model. There we estimated the (discrete) 
proportionate growth rate by using two data points, namely, the 
population sizes in 1790 and 1890. 
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This has the advantage of relative simplicity. but a more accurate answer 

can be expected if all the available data are taken into account. We shall 

now do so in the context of the continuous exponential model, to estimate 
a value for the proportionate growth rate K. 


Year 1790 1800 1810 1820 1830 1840 1850 1860 1870 1880 1890 The population is given in 
Population] 3.9 5.3 7.2 9.6 129 17.1 232 314 39.8 50.2 62.9 Millions, to the nearest 
100 000. 
We start from the equation 
P= Pre®, (3.7) 
and take the logarithm of each side, This gives 
In P = In(Poe**) ; Recall that In(e") = x, and 
= In Py + In(e**) = In Py + Kt. that In(zy) = Ina + Iny. 


Now compare this equation with the general form of the linear function 
y = mea + ec, as follows: Here m means the constant 
slope of a linear function, and 


= mer c 
y f i ; not the mass of a sample of 
z radioactive substance. 
InP = Kt + nP 
This comparison shows that, provided the exponential model is Similar comparisons were 


appropriate, we have a linear relationship between the variables In P and ¢, carried out for the radioactive 
Hence a plot of In P against ¢ should result in the data points lying on or decay model, and log-linear 
close to a straight line which crosses the (In P)-axis at In Py and has plots were drawn from data. 
slope K. Such a plot of In P against ¢ is called a log-linear plot. In that case, ¢ = 0. 


Example 3.1 US population (1790-1890) 


(a) By using a log-linear plot for the data in the table above, show that 
the exponential model is appropriate for the US population for the 
years 1790-1890. Measure time ¢ in years, from ¢ = 0 at the 1790 
census date. 


(b) Use your plot to estimate values for the proportionate growth rate 
and the initial population size Pp in the exponential model. 


(c) Use this exponential model to predict the population of the USA 
in 1950. 
Solution 


(a) The table below gives the values of the time ¢, the population size P 
and In P, to three significant figures. 


t 0 1 862006300640) «6500 60 70 80 90s 100 


iP: 39 53 7.2 96 12.9 17.1 23.2 31.4 39.8 50.2 62.9 
InP | 1.36 1.67 1.97 2.26 256 284 3.14 345 3.68 3.92 4.14 


The corresponding log-linear plot is shown in Figure 3.5. 


raight line was drawn 


The value of the parameter Po 
for the model differs from the 
actual value of the population 
at time t = 0 because the 
model value is chosen by a 
process of line-fitting that 
involyes all the given data for 
the population, 
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(b 


(c 


: 
9 20 40 60 80 100 ¢ 


Figure 3.5 Plot of In P against t for the data in the table 


The data points are fitted reasonably well by the straight line drawn in 
the figure. This confirms that the exponential model is appropriate to 
describe the US population for the years 1790-1890, 


The intercept of the line on the (In P)-axis is 1.40; this is the value of 
In P, for the model. Since In Py = 1.40, we have 
P, =e! = 4.06 (to 3 s.f.). 
As well as the point (0.1.40), another point on the straight line is 
(100, 4.23). The slope of the straight line is therefore 
4.23 — 1.40 
100 —0 


So our estimate for the proportionate growth rate K is 0.0283 year~'. 
and for the initial population size P, the estimate is 4.06 (million). 


= 0.0283. 


Using the estimates obtained for K and P), the exponential model for 
the US population (1790-1890) is 


P = 4.060928, 


The year 1950 corresponds to t = 160, so our estimate for the 
population in that year is 


*P(160) = 4.06€"°?83*!© — 376 (to 3 s.f.). 
So our estimate for the population of the USA in 1950 is 376 million. 


Comment 


The actual US population in 1950 was only 150.7 million! This shows that 
the exponential model used here cannot be expected to give good results 
much beyond the period 1790-1890 for which the data were given. 


In the following activity, you can carry out a similar exercise for another 
population. 
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Activity 3.7 Population of India (1950-2000) 


The table below shows data for the population of India (in millions, to the 
nearest million) for the years 1950-2000. 


Year 1950 1960 1970 1980 1990 2000 
Population | 358 442 555 689 851 1014 


(a) By using a log-linear plot for the data in this table, show that the 
exponential model is appropriate for the Indian population for the 
years 1950-2000. Measure time ¢ in years, from t = 0 at the 1950 
census date. 


(b) Use your plot to estimate values for the proportionate growth rate K 
and the initial population size Py in the exponential model. 


(c) Use this exponential model to predict the population of India in 2020, 


Solutions are given on page 48. 


Mathematical note 


For each of the radioactive decay and population models, we restricted the 
dependent variable (m and P, respectively) to be positive. In terms of the 
independent variable x and dependent variable y, we have derived the 
result that dy/da = Ky (y > 0) has general solution y = Ae**, where A is 
a positive but otherwise arbitrary constant. The restriction on y arose 
from modelling considerations, though it was also relied on in the niathide ‘The method used the fact 
of solution. that 


If the restriction on y is removed, then the expression above for the general if day =Iny+c, 
solution remains correct, but with no restriction on A; that is, y 
a : for y > 0. An alternat ive 
—=Ky has the general solution y = Ae*?, (3.9) approach is required to deal 
dx with the cases y < 0, 
where A is an arbitrary constant (which may be positive, negative or zero). 
This solution can be checked, as usual, by verifying that the differential 
equation is valid when y = Ae? is substituted into it. 


Comparison with the discrete exponential model The remainder of this 
This subsection concludes by considering how the continuous exponential subsection will'not be 
assessed. 


model for population change compares with the corresponding discrete 
model developed in Chapter B1. The discrete model was given by the 
recurrence relation See Chapter B1, 


Pun =(14-r)Be Gr S0meeee), SE les 


where P,, is the population size at n years after a specified initial time and 
r is the (discrete) proportionate growth rate (giving the factor by which 
the population increases year on year). This recurrence relation has the 
closed-form solution 


P,=(1+r)"Po (n=0,1,2,...). 


From either the recurrence relation or its solution, we have 


Meith, (3.10) 
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For example, if K = 0.03 then 
r ~ 0.0305, and if K = 0.2 
then r ~ 0,22. 


The discrete model does, to 
the extent that it focuses on a 
‘snapshot’ at annual intervals. 


For example, the continuous 
logistic model, which is the 
analogue of the diserete 
model in Chapter B1, 
Sections 3-4, is given by the 
differential equation 


where K and FE are constants. 


36 


CHAPTER C3 DIFFERENTIAL EQUATIONS AND MODELLING 


This can be compared with the corresponding formula for the continuous 
model, which is 

ze rilh 

P(t) 

Now P,, and P,,.; are population sizes taken one year apart, at t =n and 
t=n-+1. Hence we can match equations (3.8) and (3.10), for integer 
values of t, by putting t = n and T = 1 into equation (3.8). We see then 
that the two models will give the same values for P(n) = P,,, provided that 


KT 


(3.8) 


eX=1+r. 


This equation relates the discrete proportionate growth rate r (measured 
over one year), with the continuous proportionate growth rate K 
(measured instantaneously). For a given population, these two 
proportionate growth rates are not quite the same. However, if the 
population size varies slowly, then the two values will be close to each 
other. 


As was pointed out, a continuous model would be inappropriate for a 
population whose size is small. Continuous models have the apparent 
advantage over discrete models that they predict the population at all 
times, and not just annually. However, for many species the birth rate is 
highest in the spring whereas the death rate is highest in the winter. The 
continuous exponential model does not take account of this seasonal 
variation. 

Furthermore, as discussed in Chapter B1, for very large populations the 
environment becomes less hospitable. Birth rates then tend to decline and 
mortality rates to rise, so in practice the proportionate growth rate 
decreases with the size of the population. There are continuous models, 
based on differential equations, which take account both of seasonal 
variations and of the decrease in proportionate growth rate with population 
size. However, there is not the space to study them in this course. 


Summary of Section 3 


This section has introduced: 

© the (continuous) exponential models for radioactive decay and for 
population change; 

© a method for finding parameter values in these models from data, 
using a log-linear plot: 

© the form of differential equation dy/dx = Ky (occurring in each of the 
above models), with general solution y = Ae. 
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Exercises for Section 3 


Exercise 3.1 

The interest on bank deposit accounts is usually compounded annually, Compounding interest 
six-monthly or quarterly. Suppose that, in order to attract customers, a n times per year, at a rate 
certain bank introduces a type of account on which interest is compounded” per year, means that an 
continuously at the rate of 5% per year. initial credit balance co rises 


to 
(a) Use the fact that the proportionate rate of increase of any credit 


pyr 
balance is ;3; = 0.05 to write down a differential equation for the 0 (1 + 4) 
credit balance ¢ (in £) of such an account at time ¢ (in years), during 


: 5 . rH fer year. 
any time interval in which no deposits or withdrawals are made. nk al 
(b) Write down the general solution of this differential equation. 


(c) A customer places £100 in such an account and makes no further 
deposits or withdrawals in the subsequent year. Find the 
corresponding particular solution of the differential equation, and 
hence calculate the credit balance of the account after one year. 


(d) Compare the interest received by this customer after one year with 
that which would be obtained at a rate of 5% per year compounded 


(i) annually; (ii) six-monthly; (iii) quarterly. 
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4 Differential equations with the computer 


Ca) 


For each point (x,y) in the 
(two-dimensional) domain 
of f, a unique real value 
f(x,y) is defined. 
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To study this section you will need access to your computer, together with 
Computer Book C. 


In Section 2, you saw how to solve equations of the form dy/dxr = f(x)g(y) 
using the method of separation of variables. There are many first-order 
differential equations which cannot be expressed in this form. For example, 
dy 
Fe =ar+y (4.1) 
is one such equation. More generally, we can consider differential equations 
of the form 


dy 

ae 7 f(my)- (4.2) 
Here f is a known function of two variables, «x and y. In some cases such 
differential equations can be solved to obtain a formula (either implicit or 
explicit) for the solution, Such a solution is said be analytical. On the 
other hand, it is sometimes not possible to find a formula which describes 
a solution. In such cases, the behaviour of functions that satisfy a 
differential equation can be deduced instead by numerical means, using a 
computer. You will see how this is done here. 


We start by considering the graphical information which may be extracted 
from a first-order differential equation. A solution of equation (4.2) will 
have a graph which is a curve in the (x, y)-plane. For each point (a, y) that 
the curve passes through, the gradient of the curve has the corresponding 
value f(a, y). 


However, the differential equation specifies such gradient values 
throughout the domain S, in the (7. y)-plane, of the function f, and not 
just on any particular curve that passes through S. Hence the behaviour of 
the differential equation (as opposed to any individual solution of it) can 
be represented graphically if we have a means of visualising gradients at all 
points of the (x, y)-plane. 


The association of a gradient with each point in the domain S is called a 
direction field for the differential equation. The gradient at each point 
can be represented satisfactorily by a short line segment that has the 
appropriate gradient and is centred on the point concerned. It is not a 
practical proposition to draw such line segments at all points of S, but we 
can do so for a representative grid of points. If the spacing between grid 
points is not too large, then this will provide a good picture of how the 
differential equation behaves. 


The direction field for the differential equation (4.1) is represented in 
Figure 4.1. It can be seen that the gradients shown do correspond to the 
numerical values of f(x,y) = 2+ y. For example, the gradient is zero at 
the origin, —1 at (—2.1) and 2 at (1,1). 
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Figure 4.1 Direction field for dy/da = x+y 


The graph of any solution to the differential equation must have a gradient 
at each point which matches that provided by the direction field. You can 
probably see that such graphs could readily be sketched on top of the 
direction field diagram, ‘following the flow’, to give a rough idea of how 
solutions behave. There are also more accurate ways of constructing such 
graphs, as you will see shortly. 


Refer to Computer Book C for the work in this section. 


You have seen how solutions to a differential equation of the type (4.2) can 
be obtained numerically. In fact, each such solution corresponds to a 
particular initial-value problem, 

dy 

ae f(x,y). 
This can be solved using Euler’s method, which is described by the pair 
of recurrence relations 


y(20) = Yo- 


Yne1 =Un+hf(Tn.tn) (n=0,1,2,...). 


In these equations, h is the step size between the successive values of x at 
which solution estimates are calculated. Each calculated value y,, is an 
estimate of the corresponding ‘true solution’ y at « = x,,; that is, y, is an 
estimate of y(x,,). The sequence of estimates depends on the choice of 
both the step size h and the overall number of steps N. Decreasing h, 
while increasing N to cover the same range of x-values, leads to 
progressively improved estimates for the solution values, and with a small 
enough step size, any desired level of accuracy can be achieved. 


Try = Ty +h, 


Summary of Section 4 


This section showed how first-order differential equations can be visualised 
graphically, using direction fields, and how solutions may be obtained 
either numerically or graphically, using Euler’s method. 


The Swiss mathematician 
Leonhard Euler (1707-83) 
was hugely prolific. He 
devised the computational 
method that now bears his 
name in order to calculate the 
orbit of the Moon. It was 
used by the German 
cartographer and astronomer 
(Johann) Tobias Mayer 
(1723-1762) to construct 
Junar tables. In 1765, Mayer's 
widow received £3000 from 
the British Board of 
Longitude for the 
contribution which his tables 
made to the problem of 
determining longitude, while 
Buler received £300 for his 
theoretical contribution to the 
work. 
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5 A look ahead 


This section will not be 
assessed, 


You can check this solution 
by substitution, 


This is called the simple 
harmonic motion equation, 
and provides the simplest 
model of an oscillating 
system. Again, you can check 
the solution by substitution. 


For example, = might be the 
temperature at the point 
(x,y) of a flat metal plate. 


For example, the variation in 
temperature O(a, t) at 
distance « along a thin rod at 
time t is modelled by the heat 
equation, 
00 
Ou? ot 
where k is a constant, 
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Differential equations form a vast field of study, and only the very 
beginnings have been described in this chapter. Their importance lies in 
the frequency with which they crop up as central relationships within 
mathematical models. You will see much more of differential equations if 
you undertake further studies in mathematics. 


Equations of the form dy/dx = f(a) and dy/dx = f(x)g(y) do not exhaust 
the types of first-order differential equations for which analytical solutions 
can be found. For example, the differential equation 


ey (4.1) 


dx 
whose direction field you studied with the computer in Section 4, has the 
general solution 


y= Ae* —2z-1, 


where A is an arbitrary constant. Higher-order equations can also be 
solved analytically. For example, the second-order differential equation 


where w is a non-zero constant, has the general solution 
y = Acos(wt) + Bsin(wt), 


where A and B are arbitrary constants. Note the presence of two arbitrary 
constants here, for the general solution of a second-order equation. To 
obtain a particular solution in this case, two initial conditions are required, 
perhaps corresponding to a known position and velocity at a given time. 
This is typical behaviour; the number of arbitrary constants that feature in 
the general solution is the same as the order of the differential equation. 


Where analytical solutions cannot be found, numerical methods may again 
be employed, for any order of differential equation. Euler’s method is just 
the first in a long line of methods which can be used for this purpose. 


We have so far considered equations that involve derivatives of a function 
of one independent variable, say y = f(x). It is also possible to develop 
calculus for functions of two or more independent variables, say z = f(a, y). 
Here there are two types of first-order derivative of z which may be taken: 
with respect to x, denoted 8z/Ax, or with respect to y, denoted 02/Oy. 
These are called partial derivatives, in contrast to the ordinary derivatives 
considered in this block. Equations relating such derivatives, which again 
crop up within significant mathematical models, are called partial 
differential equations. As with ordinary differential equations, these may 
be of first or higher order. While some can be solved analytically, it is 
often the case that numerical methods are required to find their solutions. 


As a final remark, it is possible to combine these ideas with those of 
vectors, which you saw in Chapter B3, and hence to arrive at the topic of 
vector calculus. Both ordinary and partial differential equations feature 
here. 


Summary of Chapter C3 


In this chapter you were introduced to first-order differential equations. 
You saw how to solve equations of the form dy/dx = f(x) by direct 
integration, and equations of the form dy/da = f(:r)g(y) by separation of 
variables. The more general type of equation dy/dx = f(x,y) can be 
solved numerically. 


The continuous exponential model, embodied in the differential equation 
dy/dx = Ky (where K is a constant), has numerous applications. You saw 
it applied to radioactive decay and population change. If data are available 
in such a case, then values may be found for the parameters (such as K’) 
from a log-linear plot. 


Learning outcomes 


You have been working towards the following learning outcomes. 


Terms to know and use 


Differential equation. order, first-order, solution, general solution, 
particular solution, initial condition, initial values, initial-value problem, 
direct integration, implicit solution, explicit solution, implicit 
differentiation, separable, separation of variables, exponential 
growth/decay, decay constant, half-life, proportionate growth rate, 
doubling time, discrete/continuous exponential model, log-linear plot, 
direction field, Euler’s method, step size. 


Notation to know and use 
y(a) = b (to express an initial condition on the variable y). 


Mathematical skills 


© Be able to recognise differential equations which can be solved by 
direct integration, and to solve suitable equations by this method. 


© Be able to recognise differential equations which can be solved by 
separation of variables, and to solve suitable equations by this method. 


© Knowing the general solution of a differential equation, be able to use 
an initial condition to find a particular solution. 


© Know the importance of checking a solution of a differential equation 
by substituting it back into the equation. 


© Be able to write down the general solution of a differential equation of 
the form dy/dx = Ky, where K is a constant. 


© Given data for the variables in such a differential equation, use a 
log-linear plot to estimate values for the parameters (such as KX). 


© In an exponential decay process, be able to find the half-life from the 
decay constant, or vice versa. 


Modelling skills 


© Interpret predictions from an exponential model in real-world terms. 
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CHAPTER C3 DIFFERENTIAL EQUATIONS AND MODELLING 


Summary of Block C 


This block has introduced calculus, which is the cornerstone of 
mathematical studies at this level. Facility with calculus comes at the 
price of much exposure to it and practice with it, so do not be perturbed if 
you feel that you have yet to grasp fully some of the important points. 


Chapter C1 dealt with the topic of differentiation, starting with an 
informal explanation in terms of rates of change and the gradients of 
graphs, then providing a formal definition involving a limit. The 
consequences of this definition included a table of standard derivatives and 
various rules for differentiating sums and constant multiples of functions. 
products of functions, quotients of functions and composite functions. The 
alternative Leibniz notation was introduced. Differentiation was applied to 
the process of sketching graphs and to optimisation. 


Chapter C2 started from the basis that integration is both the reverse 
process to differentiation and a process in its own right. A table of 
standard indefinite integrals was drawn up, corresponding to that for 
derivatives, and rules were deduced for integrating sums and constant 
multiples of functions. The role of basic algebraic manipulation was 
stressed, as a means of modifying the expression for a function until it is in 
a form where its integral can be recognised, from the table or otherwise. 
Integration was applied to analyse one-dimensional motion, and in 
particular motion with constant acceleration, for which specific formulas 
were derived connecting time, position and velocity. It was shown that 
definite integration provides a method for finding areas under graphs. 


Chapter C3 involved in significant part revision and reinforcement of the 
material in the previous two chapters. Integration is central to the two 
strategies put forward for solving differential equations, and differentiation 
is needed in order to check solutions. (Also, in the shape of implicit 
differentiation with the Chain Rule, it permitted an explanation of why 
the separation of variables method works.) The principal modelling 
application was to exponential growth and decay, though other contexts 
were introduced also. 


This block has been but a brief initial tour of the possibilities provided by 
calculus, and you can expect to come across it again on a frequent basis if 
you continue with mathematical studies. This will provide reinforcement 
to your understanding of the basic ideas and results, as well as further 
appreciation of its use in modelling. 


Solutions to Activities 


Solution 1.1 


(a) The differential equation is dy/dx = 5/(2Vz). 
The derivative of the given function, 
y= 5 Ja = 52'/?, is 
gos 
2/r° 
which is the same expression as the right-hand 
side of the differential equation, as required. 


452") efrtee 


(b) The differential equation is 
dy/da = 3e** cos(e**). In order to differentiate 
the given function, y = sin(e**), we apply the 
Composite (Chain) Rule. The function can be 
written as 


y=sinu, where u=e*, 


for which 
dy _ du _ aa: 
iy ot and Ps =3e™. 
Hence we have 
dy _ dy du 
dx dudxr 


= (cos u)(3e**) 

= 3e* cos(e**), 
which is the same expression as the right-hand 
side of the differential equation, as required, 


(c) The differential equation is 
dy/da = y(1 — 2tan(2:r)). In order to 
differentiate the given function, y 


0s (2), 


we apply the Product Rule. The function can be 


written as y = uv, where 
u=e" and v=cos(22), 
for which 
du dv 


—=e — = —2sii 
= and Ee 2sin(2r). 
Hence we have 


= (e*)(cos(2a:)) + (e*)(—2sin(2x)) 

= e*(cos(2r) — 2sin(2r)). 
Now we substitute y = e* cos(2r) into the 
right-hand side of the differential equation, to 
obtain 

y(1 — 2tan(2r)) = e* cos(2a)(1 — 2tan(2x)) 
= e"(cos(2x) — 2sin(2z)). 

This is the same expression as was obtained 
above for the derivative, as required. 


(d) The differential equation is dy/dx = —3(y +1). 
In order to differentiate the given function, 
y=(1—2x)/(1+~2), we apply the Quotient Rule. 
The function can be written as y = u/v, where 


u=1-—2 and 1+z, 
for which 

du dv 

hd -1 and a a 


Hence we have 


dy 1 du dv 

ay 1 (tit) 
— @+2)(-1) - (1—2)() 

a (+2)? 

2 


~~ Faye 

Now we substitute y = (1 —x)/(1 +.) into the 
right-hand side of the differential equation, to 
obtain 


This is the same expression as was obtained 
above for the derivative, as required. 

Solution 1.2 

Substitute the values «= 1 and y = 7 into the 

expression y = 5% +c. This gives 
T=5VI+e=5+e, 

so ¢= 2. The required particular solution is therefore 
y=5¥r+2. 


Solution 1.3 
The direct integration method can be applied to a 
differential equation of the form 
d(dependent variable) A 
iadependentresabiey— f (independent variable). 


Equations (a), (c) and (e) are of this form, so direct 
integration can be applied to these. Equations (b), 
(d) and (f) are not of the appropriate form. 


43 


CHAPTER C3 DIFFERENTIAL EQUATIONS AND MODELLING 


Solution 1.4 
(a) The general solution of dy/da = e2* is 


[eu 


aie 
= jet +e, 


where ¢ is an arbitrary constant. 


(Check: Differentiating y = 3e”* + ¢ gives 
dy/dx = e?*, as required.) 

(b) The initial condition is y = 2 when x = 0. 
Putting « = 0 and y = 2 into the general 
solution, we obtain 


ep rea? | P, 
Q=ge +ce=5+e. 


Hence ¢ = 3, and the required particular 
solution is 


y=fe*+3. 


Solution 1.5 


The general solution of dx/dt = sint is 


= [sineae 


=—cost+e, 
where ¢ is an arbitrary constant. 


(Check: Differentiating « = —cost + ¢ gives 
dx/dt = sint, as required.) 


The initial condition is 2 = 1 when t = 0. Putting 
t=0 and x = 1 into the general solution, we obtain 


1=-cos0+c=-l+e. 


Hence c = 2, and the solution of the initial-value 
problem is 


a= 2—cost. 


Solution 1.6 
(a) The general solution of dm/dt = —8t is 


m= [csoa 


=-4t? +¢, 
where c is an arbitrary constant. (This general 
solution can be checked by differentiation.) 


(b) The initial condition is m= 200 when ¢ = 0. 
Putting ¢ = 0 and m = 200 into the general 
solution, we obtain 


200=0+e. 


Hence ¢ = 200, and the required particular 
solution is 


m = 200 —4¢?. 


(c) When m = 100, the corresponding value of t is 
given by 
100 = 200 —4#?: that is, t? = 25. 


Hence ¢ = 5 when m = 100. The rockets fuel is 
exhausted after 5 seconds. 


Solution 2.1 


(a) The given equation is of the form H(y) = F(x), 
with H(y) =y+siny and F(x) = 2+ e* —1, It 
follows from result (2.3) that 


meee: Majed | oe 
al nw) Ge = gee te Ds 
which leads to 


(1 +cosy) =1+e*. 


viding through by 1 + cosy (which is non-zero 
for —m < y < 7) gives 


ay Dre 
dx 1+cosy’ 
as required. 


(b) At (0,0), the slope of the curve is 
dy l+e? 2 


= ae 


dz 1+cos0 2 


Solution 2.2 


The given equation is of the form H(y) = F(x), with 
H(y) = tany and F(z) = 2. It follows from 
result (2.3) that 
d dy a 
ga Pri 
which leads to 


dy 
ey =1. 
sec” y 1 


ing through by sec? y (which is non-zero for 
—}7 <y < $7), and using the trigonometric identity 
sec? y = 1 + tan? y, we have 

dy 1 

dx 1+tan?y" 
However, tan y = -r, and so 


dy 1 
dr 1+2?° 
as required. 


Solution 2.3 


(a) 


(b) 


The differential equation is (1/y*)dy/dx = —1. 
Comparing with result (2.4), we have 

h(y) = 1/y? = y~? and f(x) =—1. Hence the 
general solution of the differential equation is 

given by 


[vra= [nen 


that is, 
--=-1+e 


where ¢ is an arbitrary constant. 


The explicit solution is obtained by making y the 
subject of the above equation. In this case, that 
is achieved by taking the reciprocal of each side 
and multiplying through by —1. The outcome is 


1 


y= 


where ¢ is an arbitrary constant. (Since y > 0, 
we must have « —c¢ > 0. There is some 
discussion of such restrictions on variables later 
in the section.) 


This answer could equally well be written as 
1 
rte’ 


where ¢ is an arbitrary constant (in which case, 
we must have x + ¢ > 0). 


y= 


Solution 2.4 


Equations (a), (b), (c) and (e) are of the appropriate 
form, but equations (d) and (f) are not. 


(a) 
(b) 
(c) 
(e) 


f(z) =2, gly) = cosy 
f(t) =teost, g(x) =1 
F(t)=1,  g(p) = peosp 
f(y)=1+y. glv)=v 


These are the most likely choices for f and g, but 
they are not the only correct possibilities. For 
equation (a), for example, you could have chosen 
f(a) = 2x and g(y) = } cosy. 


Equation (b) can also be solved using direct 
integration; this is the most direct method of 
solution, 


SOLUTIONS TO ACTIVITIES 


Solution 2.5 


(a) 


(b) 


The differential equation dy/dxz = —a/y has a 
right-hand side of the form f(r)g(y), where 
f(x) =—wx and g(y) = 1/y. Dividing both sides 
of the differential equation by g(y) means 
multiplying both sides by y. This gives 
d 

ve =z. 
(Other choices for the functions f and g are 
possible, but all choices lead to the same general 
solution.) 


From this we obtain, using result (2.4), 


f[vav= [oae. 


(The solution thus far can be shortened by using 
the shortcut described in the text before 
Activity 2.4 on page 20. In that case, the 
solution to this point would simply read: 
Separating the variables in dy/dz = —x/y, we 
obtain 


[ vtv= [(-oae) 


Carrying out the two integrations, we have 
3y? =-}2" +c, 


where ¢ is an arbitrary constant. This implicit 
form of the general solution can be manipulated 
to make y the subject, giving 


y= V2e—2?, 
or equivalently 

y= vb—22, 
where b = 2c is another arbitrary constant. We 
take the positive square root for y in order to 


satisfy the restriction y > 0 given in the 
question. 


(This expression for the general solution could 
be checked by substituting it into the differential 
equation.) 

The initial condition is y= 1 when x = 0. 
Putting 2 = 0 and y = 1 into the general 
solution, we have 


inane 


Hence 6 = 1, and the required particular solution 


Is 
y= Vi1-2. 


(For this expression to be well-defined, the range 
of values of x: must be restricted. We shall 
return to this problem later in the section.) 
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Solution 2.6 


(a) 


(b 


) The shortcut is adopted here. The differential 
equation is dir/dt = t sec x, so we can choose 
f(t) =t and g(x) =secx. Separating the 
variables in the differential equation, we obtain 


|meeee [eae 
sec 


that is, 


[ cosivae = [ta 


After integration, we have 


sine =}? +e, 


where ¢ is aa arbitrary constant. Since 
- 3m <a< hr, making the dependent variable x 
the subject ii the equation, we obtain 


= aresin($t? + ¢) 


for the general solution in explicit form. 


The initial condition is « = 0 when t = 0. 
Putting ¢ = 0 and x = 0 into the general 
solution gives 


0 = aresin(0 + ¢) = aresine, 


Hence c¢ = 0, and the required particular solution 
is 


x = aresin($#*). 


Solution 2.7 


(a 


(b 


) Separation of variables in the differential 
equation gives 


fa hee 


Hence we have 


k 
= ste 


or equivalently (since v > 0) 


k 
v= 1/2 ( ne °). 
r 
where ¢ is an arbitrary constant. This is the 


general solution of the differential equation. 


) The initial condition is v(R) = va. Putting 
r= R and v = vp into the general solution, we 
have 


viata) 


Hence c = $v} — k/R, and the required 
particular solution is 


1 
12, -<—— 
o=yeirai (lag) 


46 


(c) As r becomes large. 1/r — 0, so the velocity 


v vz +2k 


(:-#) 


tends towards the limiting value 


If vp > \/2k/R then vj > 2k/R, and so 


2k 
V°R- x re 


Therefore, if uz > \/2k/R, the probe continues 
to travel away from the Earth indefinitely, since 
its velocity is never zero. 


On the other hand, if up < \/2k/R then 

v}, < 2k/R. and so no real value for the velocity 
is obtained at sufficiently large distances from 
the Earth. In this case the probe does not travel 
that far away; it attains a maximum distance 
from the Earth and then falls back. 


(An alternative approach to this question would 
be to examine where the velocity of the probe 
becomes zero. This occurs where 


ra(-h)” 
cake 
which gives a positive value for r only if 


up < \/2k/R. Hence the velocity of the probe is 
never zero if up > \/2k/R.) 


The critical launch velocity, 
called the Earth’s escape 
significance is that, according to the model, a 
probe launched with this or a greater velocity 
will escape completely from the Earth, whereas 
a probe launched with a lesser velocity will 
eventually fall back towards the Earth. (The 
value of the Earth’s escape velocity is 
approximately 1.12 x 10'ms*!, that is, about 
40000 kph or 25000 mph.) 


Solution 3.1 
(a) The differential equation is 


dm/dt = —km (m > 0). Separating the 
variables, we obtain 


ee 
m 


Since m > 0, this gives 
Inm = —-kt+e, 


where ¢ is an arbitrary constant, This is the 
general solution in implicit form. 


SOLUTIONS TO ACTIVITIES 


Taking the exponential of each side, we obtain Solution 3.4 
mae (a) The completed table of data is below. 
aie = Ae, 
t (min) 10 20 30 40 50 
where A = e* is a positive but otherwise - 
arbitrary constant. This is the general solution m/mo O73 0.56 042 0.29 0.24 
in explicit form. In(m/mo) | —0.31 —0.58 —0.87 -1.24 -1.43 
(b) The initial condition is m= mo when t = 0. ‘ . 
Putting t = 0 and m = mo into the general (b) The graph of In(rn/mo) against t for the data in 
solawonswelhave part (a) is shown in Figure S.1. This figure also 
7 shows a straight line (chosen by eye) through 
mo = Ae” = A. the origin which fits the data points well. (This 
Hence A = mo, and the required particular good fit confirms that the radioactive decay of 
splutiea ts uranium-239 can be modelled satisfactorily as 
= an exponential decay process.) 
m= moe*. 
tal /me} 10 20 30 40 50 
Solution 3.2 


Using equation (3.5), the decay constant k 
corresponding to a half-life 7’ = 2.62 hours is 


In2_ In2 =i . 
7 *Te7 0.265 hour™! (to 3 s.f.). 
Hence, from equation (3.4), the proportion of the 
original amount remaining after six hours is 
m(t+6) _ -o.265%6 _ 

may oe = 0.204 (to 3 s.f.). 
So about 20% of the original amount of silicon-31 
remains after six hours. 


k 


Solution 3.3 
Using equation (3.5). the decay constant k 
corresponding to a half-life T = 5570 years is Figure S.1 
(c) The straight line passes through the origin and 
At the present time, 85% of the original amount mo the point (50, —1.48), 80 its slope is 
of carbon-14 remains, so we have 8 = ~0.0296. 
at, Sekt 
8b mo en Hence the decay constant for uranium-239 is 
a in72 ion (¢ 4 
where ¢ (in years) is the time since the animal died. be a ese on equation (3.5), the 
Taking the logarithm of each side, we obtain Se 
In2 In2 A 
—kt = In(0.85); T = = = pppgg = 23-4 min (to 3 sf), 
that is, 
_ 1n(0.85) 
xe k 
__1n(0.85) 


"= 1310 years (to 3 s.f.). 


So the animal died approximately 1300 years ago. 
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Solution 3.5 


(a) The initial population size is Py = 5 x 10°, and 
the proportionate growth rate is 
K = 0.03 year~'. Hence, from equation (3.7), 
the population size P after time ¢ years is 


P= PreXt = 5 x 10%, 
After 10 years the population size is 
P=5 x 10%" = 6.75 x 10° (to 3s.f.), 
whereas after 100 years, the population size is 
P=5 x 10%% = 1.00 x 10° (to 3 sf). 


So, to three significant figures, the population 
size after 10 years is 6.75 million, and after 100 
years it is 100 million. 

(b) The population size reaches 50 x 10° at the time 
t years given by 


50 x 10° =5 x 10%"; that is, e°° = 10. 


Taking the logarithm of each side, we obtain 
0,03 = In 10; that is, 


t = ln 10 = 76.8 (to 3 s.f.). 


Hence it takes almost 77 years for the 
population size to reach 50 million, 


Solution 3.6 


(a) From equation (3.8), the population will double 
in a time T given by 
Pe+T) _ kT 
Py 
Taking the logarithm of each side, we obtain 
KT = \n2, so the doubling time is 
In2 
T= K 
(b) For the population described in Activity 3.5, the 
proportionate growth rate is K = 0.03 year™!. 
Hence the doubling time is 
nz 
0.03 


2= 


T= = 23.1 years (to 3 s.f.). 


Solution 3.7 


(a) The table below gives the values ofthe time ¢ 
(in years since 1950), the population size P (in 
millions) and In P, to three significant figures. 


t 0 10 8620 6300640050 
P 358 442 555 689 851 1014 
InP} 5.88 6.09 6.32 6.54 6.75 6.92 


The corresponding log-linear plot is shown in 
Figure S.2. 
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(b) 


(c) 


Figure S.2 


‘The data points are fitted reasonably well by the 
straight line drawn. This confirms that the 
exponential model is appropriate to describe the 
Indian population for the years 1950-2000. 

The intercept on the (In P)-axis is 5.89; this is 
the value of In Po for the model. Since 

In Po = 5.89, we have 


Py = e*** = 361 (to 3 sf.) 


As well as the point (0.5.89), another point on 
the straight line is (50,6.94), The slope of the 
straight line is therefore 

6.94 — 5,89 

50-0 

So our estimate for the proportionate growth 
rate K is 0.021 year~!, and for the initial 
population Po the estimate is 361 (million). 


= 0.021. 


Using the estimates obtained for K and Po, the 
exponential model for the Indian population 
(1950-2000) is 


P= 361", 


The year 2020 corresponds to t = 70, so our 
estimate for the population in that year is 


P(70) = 361¢°°7!*79 = 1570 (to 3 s.f.). 


So our estimate for the population of India in 
2020 is 1570 million, that is, 1.57 billion. 


Solutions to Exercises 


Solution 1.1 


The differential equation is dy/dx 
derivative of the given function, y = 


£(2)~-dere+ 


Now we substitute y = —1/ into the right-hand side 
of the differential equation, to obtain 


2 2 (peas 
u -( zt)” x" 


This is the same expression as was obtained above 
for the derivative. Hence y= —1/z (x > 0) isa 
solution of the given differential equation. 


Solution 1.2 
(a) The general solution of dy/dr = = 


rf 


=}: 


al? js 


? dar 


3/2 +c, 
where ¢ is an arbitrary constant. 


(b) The initial condition is y= 5 when x = 4. 
Putting « = 4 and y = 5 into the general 
solution, we obtain 


cons 


xPPte= Bie 


Hence ¢ = 
solution is 


. and the required particular 


y= 9c? — 4. 
(c) When x = 25, we have 
x 259/221 
x 125-4 
= 4(250 — 1) = 83. 


y= 


esine coins 


Solution 1.3 


(a) The general solution of du/dx = cos(2x) is 
c= costa) ae 


= fsin(2x) +c, 


where ¢ is an arbitrary constant. Using the 
initial condition, u = —2 when x = dn, we 
obtain 


—2= }sin(dn) +e=}+e. 


Hence c = ~$, and the solution of the 
initial-value problem is 


u= 4sin(2r) — 3. 


(b) The general solution of dx/dt = 1/t (t > 0) is 


1 
ew eae 
r fi it 
=Int+e, 


where c is an arbitrary constant. Using the 
initial condition, z = 3 when t = 1, we obtain 


3=Inl+e=c 


Hence c = 3, and the solution of the initial-value 
problem is 


a=Int+3. 


Solution 1.4 


(a) The general solution of the differential equation 
kis 


where ¢ is an arbitrary constant. 


(b) Since the initial radius of the snowball is 100cm, 
the initial condition is r = 100 when t = 0. 
Substituting this into the general solution, we 
obtain ¢ = 100. Hence the required particular 


solution is 


r = 100—kt. 


(c) The initial volume of the snowball (in cm’) is 
4x x 100%. If we denote the radius of the 
snowball after 2 days by R, then 

s 4 w 
4nR® = } x $7100", 
from which we find 
R=2-"8 x 100, 


Using the condition r = R = 2~'/3 x 100 at 
t = 2 in the particular solution of the differential 
equation, we obtain 


2-1/3 x 100 = 100 — 2k. 
It follows that 
k=50(1 — 27/4) ~ 10.315. 


(d) The snowball disappears when r = 0; that is, 
when 100 — kt = 0. This occurs when 

100 2 

e Ta2e 


= 9,695. 


It therefore takes about 9.7 days for the 
snowball to disappear. 
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Solution 2.1 


(a) The given equation is of the form H(y) = F(x). 
with H(y) =y+Iny, where y > 0, and 
F(x) = (1 — 2x)e**. It follows from result (2.3) 
that 
d dy_d ea 
= = =—((1- R 
ay” +Iny)7- a! 2r)e**), 


which (using the Product Rule on the right-hand 
side) leads to 


(: + *) a — (ay 08) + (1 —22)(3e%); 
y) dx 
that is, 
y+ dy aright 3r 
aS ge = Bae. 
Multiplying through by y/(y +1) gives 
dy _ y(1—6x)e* 
dx oy +1 
as required. 
(b) At (0,1), the slope of the curve is 
dy _A(1=O)e® 1 
de 14100 2 
Solution 2.2 
(a) The differential equation is u? du/dx = x*. 
Comparing with result (2.4), we have h(u) = u? 
and f(«) = 2*. Hence the general solution of the 
differential equation is given by 


if wdus / ade; 


where ¢ is an arbitrary constant. Solving this 
equation for u gives the explicit form of the 
general solution as 


u= (32° +5) 


where b = 3c is an arbitrary constant. 


1/3 


(b) The initial condition is v= 1 when «= 1. 
Putting « = 1 and u = 1 into the general 
solution, we obtain 

+0) 1/3 
Hence b = 2, and the required particular 
solution is 
u= (32°42). 
Solution 2.3 


(a) The differential equation is 
dy/dx = —y' (y > 0). Separating the variables, 
we obtain 


(-B)- fom 


50 


(b 


that is, 


aS 

ae =zr+e, 
where ¢ is an arbitrary constant. To make y the 
subject of this equation, first take the 
reciprocals of both sides and then divide through 
by 2, to obtain 
Fees ey 

a+) 
Since y > 0, the general solution in explicit form 
is therefore 


y 


= 
V2e+0)" 
The initial condition is y= } when x = 0. 
Putting x = 0 and y = 3 into the general 
solution, we obtain 

"ae 

27> V2 


Hence ¢ = 2, and the required particular solution 


The differential equation is dx/dt = tx (x > 0). 
Separating the variables, we obtain 


[pace [eae 
r 


that is, 


142 
Inc =P +e, 


where ¢ is an arbitrary constant. To make x the 
subject of this equation, recall that 

exp(Ina) = x. Hence, taking exponentials of 
both sides, we have 


2) 
z=el@ are 


= eve /? = Ae /?, 


where A = e*. The constant A must be positive 
(since e® > 0 for any choice of ¢) but is otherwise 
arbitrary. The general solution for 2 > 0, in 
explicit form, is therefore 


x= Act? 
where A is a positive but otherwise arbitrary 
constant. 


The initial condition is x = 2 when t = 0. 
Putting t = 0 and x = 2 into the general 
solution, we have 


2= Ae’. 
Hence A = 2. and the required particular 
solution is 


Hr 
z= et 7, 


Solution 2.4 


The volume V is decreasing, and the rate of change 
of volume dV/dt is proportional to the surface 
area A. In other words, 


dv 
=~ =—kA, 
dt c 
where k is the (positive) constant of proportionality. 


Putting V and A in terms of r. this differential 
equation becomes 
d 43 2 
a (3ar3) = —4zkr?. 
Hence, using the Chain Rule, we have 
inthe Se 
dt 
which gives 
dr _ 


Solution 2.5 


(a) The differential equation is 
dh/dt = —kh'/* (h > 0). Separating the 
variables, we obtain 


[rear = [eryar 


that is, 
2h'/? = —kt+e, 
where ¢ is an arbitrary constant. Hence the 
explicit form of the general solution is 
h='e-kt). 
(b) The initial condition is h = ho when t = 0. 


Putting t = 0 and h = ho into the general 
solution, we obtain 


Now the equation 2h'/? = —kt + ¢ above shows 
that ¢ must be positive, since we have h'/? > 0, 
k>0O and t > 0. Hence c= 2hj/?, and the 
required particular solution is 


h= (ni? — 41)”. 


(c) The tank empties to the level of the hole when 
h=0. This occurs when 


ni? — 


that is, at time 


ani? 


t= 


SOLUTIONS TO EXERCISES 


Solution 3.1 


(a) 


(b 


(c 


(d) 


The credit balance of the account increases at a 
proportionate rate 0.05, so the balance c at 
time ¢ satisfies the differential equation 

de 

—, = 0.05c. 

dt 
By comparison with result (3.9), the general 
solution is 

c= Aco t, 
where A is an arbitrary constant. (However, we 
require A > 0 for the purposes of the model.) 


If the initial balance is £100, then the initial 
condition is ¢ = 100 when t = 0. Substituting 
these values into the general solution, we obtain 
100 = Ae” = A, 
so the required particular solution is 
c= 1009, 
The credit balance after one year is therefore 
e(1) = 100e°° = 105.13, 
to the nearest penny. At the end of one year, the 
balance of the account is £105.13. 
In this part we use the formula given in the 
margin note, with co = 100. 


The account balances (in £) after one year, for 
an interest rate of 5% per year (r = 0.05) 
compounded annually (n = 1), six-monthly 
(n = 2) and quarterly (n = 4), are respectively 
as follows: 

100 x 1.05 = 105; 

100 x (1.025)? = 105.06; 

100 x (1.0125)* = 105.09. 
The interest accrued on £100 after one year at 
5% interest is therefore £5 (compounded 
annually), £5.06 (compounded six-monthly), 
£5.09 (compounded quarterly) or £5.13 
(compounded continuously). 
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